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5 : Abstract 

JL . We investigate the errors due to the use of unphysical values of light quark masses 

^ \ in lattice extractions of A functional form for the pion mass dependence of 

the quarkonium mass splittings (Am) is given as an expansion in m.,^/ (47r/,r) 
and mT^rB, where is the quarkonium Bohr radius. We find that, to lowest 
order, Am ~ A + Bml, where the scale of B is given by f^r^- To order m^ 
there are four unknown coefficients, however, utilizing multipole and operator 
product expansions, symmetry arguments eliminate one of the four unknowns. 
Using the central values for the lattice spacings which were extracted using two 
different, unphysical values for the pion mass, we find that the errors introduced 
by extrapolating to the physical regime are comparable to the errors quoted due 
to other sources. Extrapolation to physical values of the pion mass increases the 
value of as{Mz), bringing its value closer to the high energy extractions. 



1 Introduction 



Precision measurements of as are motivated not only by our need for more precise predictions 
within the confines of perturbative QCD, but also by the possibility that such a measurement 
could lead to hints of new physics. For instance, GUTS often make predictions for the value 
of as at low energies, and thus, these theories may be tested via a precise determination of as- 
Presently, there seems to be a slight discrepancy between low and high energy extractions of 
as m , and it has become crucial to have a good handle on the errors involved before we can 
determine whether or not this discrepancy is a signal for new physics. One of the low energy 
extractions, and the one with the smallest quoted error bar, is performed via a measurement 
of the splittings between quarkonia levels on the lattice. 

There have been several lattice extractions of as(M2)0 0], with the smallest quoted error 
being given in |^. The extractions are performed by determining the lattice spacing via a 
measurement of the splittings between different quarkonia energy levels, the values of which 
are known from laboratory experiments. Unfortunately, due to computational difficulties, 
the measurements are performed in a theory with unphysical values for the masses of the 
light mesons. The values of as extracted in these unphysical theories are then assumed to 
differ little from the physical values. The justification for this assumption is that as long as 
the light quark masses are small compared the typical momenta in the bound state, the level 
splittings should be insensitive to the quark masses. 

In this note we will determine the size of these errors in a systematic expansion in the pion 
mass. More specifically, we address the errors incurred by extrapolation to physical values 
of the light quark masses and the assumption of SU (3) fiavor symmetry. We use previously 
developed techniques [|] to determine an appropriate effective field theory to describe 
the interaction of light mesons with heavy QQ bound states. Once the correct effective field 
theory is found we are able to calculate the errors introduced due to extrapolation to physical 
values with no guess work on unknown functional forms. 

The remainder of this paper is structured as follows. In the next section we will review 
the lattice measurements with emphasis on the approximations used. In section three we 
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present our results of the errors. In section four we write down the chiral Lagrangian which 
describes the interaction of the quarkonia with the hght degrees of freedom. In section five 
we show that, within the confines of a multipole expansion, it is possible to map the theory 
of heavy quarks onto a chiral Lagrangian using symmetry arguments first noted in The 
final section is reserved for conclusions and the need for further computations. 

2 Present Lattice Extractions 

The accuracy of present lattice calculations are limited by the inability to properly simulate 
the light quarks. In general, it is difficult enough to include dynamical fermions in the calcu- 
lations, no less calculate using their physical masses. Thus, to measure physical parameters 
with accuracies which are competitive with experiment, it is necessary to look for observ- 
ables whose values are not strongly dependent on the dynamics of the light quarks. This is 
one of the primary reasons why heavy quark anti-quark bound states have been chosen as 
laboratories for measuring the value of a^. The level splittings in quarkonia are expected to 
be insensitive to the masses and the number of light quarks due to the fact that the Bohr 
radius is so small compared to the hadronic scale 1/Aqcd- It is therefore hoped that cal- 
culating with unphysical masses for the light mesons should be a very good approximation, 
and the computational problems which arise in the physical region of parameter space can 
be avoided. 

In this section we review the method by which the measurements are made following 
the work of ref [^. There are several sources of errors in these calculations, all of which are 
claimed to be under control. Here we will focus only on those errors induced by extrapolation 
to physical values of the light quark masses. 

The coupling is measured by first determining the lattice spacing, a, in a theory with 
some unphysical value for the pion and kaon masses. This is done by measuring the level 
splittings in quarkonia and then extracting the value of a by assuming that the value of the 
splitting in this unphysical theory is very close to physical value. Once the lattice spacing is 
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Table I: Measurements of ap{8.2 GeV). 





1S-2S 


1-S-lP 


aruq = .01, Uf = 2, am^ = 0.419(2) 


0.1793(16) 


0.1777(23) 


amq = .025, nj = 2, am^ = 0.269(1) 


0.1760(35) 


0.1735(28) 


nf = 


0.1551(11) 


0.1505(9) 



known, the coupling ap is determined by its plaquette value. 

47r 3 41 

- In 1^1,1 = —ap{ ) [1 - (1.185 + .070n/)ap] . 

o Qj 



The value of in the MS is then determined through the relation 

(2) 



1 + +Oial'f 



The measurements are performed with varying numbers (up to n/ = 2) of light quarks 
with several different, unphysical quark masses. The latest published results from the 
NRQCD collaboration are shown in table The physical values of as are then found 
by extrapolating to n/ = 3 using the fact that the inverse coupling is known to be 

almost linear for small changes in Uf. The assumption is made that the results are inde- 
pendent of the pion mass and SU{3) symmetry breaking effects are ignored. No errors are 
quoted for these parts of the procedure. 

The value of the coupling at the scale Mz is then determined by first running down to 
the mass of the charm quark using the three flavor beta function and then running back up 
to Mz, taking into account the b quark threshold. Errors due to variations in the thresholds 
were shown to be negligible. Using m^a = .01 the values found were: 

aj^{Mz) = 0.1152(24) 1^ - IP (3) 

ajjs{Mz) = 0.1154(26) 1S-2S. (4) 

The errors quoted here were found assuming a coefficient of one for the ctp term in (2). 

3 Results: Procedure for Calculating the Errors 

We will now present the results for the errors which will be derived in a later section. In 
a QCD like theory with an arbitrary pion mass, the dimensionless mass splitting can be 



written as 

1 



aAm = aA + —f 



a3' 



(4vr/) 



2 



[am. 

2 



,)Mog^+3(a™.)Mog^ + 4(a™,-)^.og^) 
(3(am,)'' + 4(amK)'' + (am,,)*) + .... (5) 



Here a is the lattice spacing, A, B, C, D are dimensionful unknowns, and the higher order 
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terms left out are suppressed by powers of ^j-- We will show in a later section that within 
the confines of a multipole expansion, the coefficient C is 1 at leading order under certain 
conditions. Furthermore, in general, it is be possible to calculate B using potential models, 
but we will leave it as an unknown for now. The parameter / is independent of m^r, and its 
value can be calculated in terms of m^r and /^r in QCD. In a given theory with light quark 
mass rriq, we know only the dimensionless combination am^r. In the two theories which have 
been investigated to date am.n- = 0.269(1), 0.419(2) for aruq = 0.01,0.025 respectively. 

Furthermore, all the simulations which have been performed with dynamical fermions 
have only included two light quark flavors. We therefore reduce the general form (^ to the 
case of SU{2) 



aAm 



Bo x9 / (amTr)^ 1 , ml \ Do(a) , 

aA2 + ^/2(am.)2 -2 + V^TTTTlog^ +T7^(«^-) 



(6) 



(a/)2 167r2 °47r/i2y Von'^a^ 

In the future, when simulations are performed with three light quarks, we will be able to 
determine the corrections due to SUi^ breaking using (|^). Notice that going to the SU{2) 
case does not reduce the number of unknowns. 

Since we have only two data points to fit, let us begin by considering the corrections 
only up to O(m^). In this case, we may solve for A and B and subsequently determine the 
error induced by extrapolating to physical masses for the pions. For the case of three fiavors 
one should choose the value (m^ + 2m|^)/3 for the "physical pion mass". Given that in the 
extraction process we find ap with two fiavors then extrapolate to three, we see no reason 
why this should remain the correct value. Since we are calculating errors here, we believe 
that the proper choice should be somewhere between the 140 and 410 MeV. 

Using the physical values ^""'^-^^ = 2.86 and ^I^M^ = 4.00 for the 15- IP and 25-15 



Table II: Measurements of a (91.2 GeV). 



(MeV) 


1S-2S 


IS-IP 


910 


0.1134 


0.1149 


640 


0.1154(26) 


0.1152(24) 


410 


0.1174 


0.1167 


140 


0.1184 


0.1174 



splittings in the Upsilon system respectively, we find for = 140(410) MeV 

a-' \phys = 2.56 (2.47) 1^ - IP (7) 

\phys = 2.53 (2.52) 1^ - 2S. (8) 

The authors of § found the value ap(8.32) = 0.178 (15-lP) and ap(8.08) = 0.179 (15-25) 
for = 640MeV, leading to, up to quartic terms in the pion mass, ap(8.73) = 0.178 
(IS* — IP) and ap(8.63) = 0.179 (15* — 25*), when the pion mass takes on its physical value. 
We have assumed here that InWi^i is independent of the pion mass given that it is a short 
distance quantity. Varying the pion mass was shown to change ap by less than 0.2% [Q. 
Table II shows our results for the dependence of a_{Mz) on the value of the pion mass. 

In calculating the size of the error here we must consider the fact that we have assumed 
that the expansion (^ is well behaved. Indeed, in calculating the values of the coefficients we 
have used the data for a theory in which the pion mass is on the order ~ 900 MeV, therefore 
these results should not be trusted at a quantitative level but should be a qualitative estimate 
of the error. To get an accurate value for the constants ^2,-82 and D2 it is important that 
at least one other lattice simulation be performed with light quark masses smaller than 
amq = 0.025. 

We may make the guess that the non-analytic piece dominates the counter-term as is 
sometimes done when working with chiral Lagrangians. Taking Airfi'^ to be the lattice spacing, 
we find that the results shift little. The net effect of the logs is to decrease the value of as{Mz) 
at the level of 0.5%, which is to be compared to the 3% effect found when just the piece 
quadratic in the pion mass is kept. 
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4 Chiral Lagrangians and Heavy Quarks 

Chiral Lagrangians for heavy-light systems have been utihzed for several years [0. For 
these systems there is only one scale in the theory besides the hadronic scale, namely the 
mass of the heavy quark. This scale is trivially removed by rescaling the mesonic field as was 
originally done in heavy quark effective field theory. One redefines the heavy meson field 
according to 

H{x) = e-*™«"-"i/"(x) (9) 
and treats this field as a classical static source labeled by its velocity, v. In so doing, spin 



symmetry as well as flavor symmetry (between D and B mesons) becomes manifest [pJ]j . 
Once the heavy quark mass has been scaled out, there is no question as to how the operators 
scale. 

In the case of heavy-heavy systems, chiral symmetry has been utilized in decays for 
quite sometime 0] [|13[ fl^ . The chiral Lagrangian for these systems, while equivalent 
to the current algebra approach utilized in the above references, organizes the expansion 
perhaps more naturally. Such Lagrangian have only been written down in the literature 



more recently [16||17|. Though the heavy-heavy system does not posses a flavor symmetry, 
the spin symmetry remains, contrary to what is claimed in [0. Spin symmetry breaking 
effects are suppressed by powers of the relative velocity of the heavy quarks instead of inverse 
powers of the mass as in the case of the heavy-light systems. The velocity scaling rules for 
the factorized heavy quark matrix elements will be the same as those derived for the NRQCD 
formalism [Q. For our calculation, the spin symmetry will not be relevant and therefore it 
will not be manifested in our phenomenological Lagrangian. 

Let us consider the chiral Lagrangian for the spin one 15* and 25* states. For the IS* state 
the lowest order Lagrangian, which is invariant under chiral symmetry is given by 



■'int 



+ c^h^'''^^'h^^'^*Tr{vdT)vdT.), (10) 
and similarly for the 25" state. Here, and throughout the rest of this paper, we follow the 



notation and normalizations in |2^. S is a unitary 3x3 matrix which contains Goldstone 



octet fields. ci^2,3 are dimensionful parameters whose scale at this point is unknown but will 
be determined in the next section, h^^ is the spin one heavy quarkonia field, labeled by its 
velocity f , and satisfies 

ivdh^;;^ = Q, v^h^;;^ = Q. (ii) 

We will ignore terms which are off diagonal in the quarkonia fields as they will give subleading 
contributions to the shift in the level splitting. It is possible to make the spin symmetry 
manifest by including the r] state along with the vector IS* state as is done for the heavy- 
light system, but for our purposes this is unnecessary. The leading order chiral symmetry 
breaking piece of the Lagrangian is 

Wsi, = c,h^''^''h^;>Tr{M{J: + St)), (12) 

and M is the diagonal quark mass matrix. As will be shown below, higher dimension terms 
in both the chirally symmetric as well as chiral symmetry breaking pieces of the Lagrangian 
will be suppressed by powers of 47r/. 

4.1 Multipole and Twist Expansions 

It was pointed out by Gottfried [|T^ that interactions of long wavelength gluons with quarko- 
nia should be well described by a multipole expansion, which yields an expansion in E/{r^^) 
where E is the external gluonic energy scale and is the Bohr radius. For decay processes 
this leads to an expansion in the relative quark velocity v. We are interested in the light 
quark mass dependence of the mass splittings, and therefore the expansion parameter be- 
comes m7r/(r^^). To implement the expansion we assume that the decay goes through a two 
step process. Which is to say, the hadronization process factorizes from the decay process. 
Thus, we split the Hamiltonian up as follows 

H = Hq + Hg + Hint, (13) 

where Hq acts on the heavy quarks only and includes the attractive and repulsive Coulomb 
potentials for the singlet and octet states respectively. Hg acts on the gluonic degrees of 
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freedom, and Hmt describes the interaction between the quarkonium and the hght degrees 
of freedom. We treat Hint as a perturbation and the eigenstates of the leading order Hamil- 
tonian are of the factorized form 

I ^) =1 0) I G). (14) 

I G) corresponds to the state of the dynamical gluons and is the ground state | 0) when the 
quarks are in a relative color singlet state. The coupling to the spin will be higher order in 
a velocity expansion. 

The calculation of the matrix element for two gluon emission then goes through much in 
the same way as the well known calculation of Rayleigh scattering in QED. Here, we shall 
simply state the result for the Euclidean space amplitude and refer the reader to for 
details. 

M = -^H,.\f.E' ^^_]_^„ f.E'\M, (15) 

The wave function of the state in which the heavy quarks are in a relative color singlet 
state, (ps, is an eigenstate of the Hamiltonian containing the singlet part of the Coulomb 
potential and, Dq acts on the gluonic Hilbert space. The factor l/{Ha — e — D^) accounts 
for adjoint-state propagation between gluon emissions, e is the eigenenergy of 0, and Ha is 
the repulsive adjoint Coulomb Hamiltonian acting on the quark piece of the Hilbert space. 
This so called "double-dipole" amplitude is the leading term in the multipole expansions. 

The fact that the chromo-electric field E'^ no longer depends upon the relative separation 
of the quarks and anti-quark allows us to write this amplitude in an OPE-like form where 
all the details of the bound state are in the Wilson coefficient Cn- For S wave states we may 
write 



M = - J2 Cr 



3 

nfco ' B 



n=2 even 

where 



(16) 



IGtt r Sk 1 



" J (27r)3 3 



2 



[Ha/ to + e/eo 



^ (17) 



n-l ' 



16n fg'^NY 
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and we have neglected spin symmetry breaking corrections. We may therefore reproduce 
the interactions of long wavelength gluons with the quarkonia via a (Minkowski) interaction 
Lagrangian 

oo _ -1 

Lint = E Cntl-''rlh^'^^h^;:^*{-E ■ (tD'r-'E) + ... (19) 

n=2 even 

where the ellipses denote higher multipoles as well as higher orders in a^. This sum looks 
much like the leading twist expansion in deep-inelastic scattering, except the corrections 
here are not power suppressed. Thus we will refer to the expansion in the number of field 
insertions (read coupling expansion) as the twist expansion^. 

If one assumes that the leading term in (|T9|) dominates, then it is possible to map this 
effective Lagrangian onto a chiral Lagrangian using symmetry arguments, as will be done in 
the next section. However, there is no a priori reason why this should be a good approx- 
imation. For decay processes, this expansion may not be useful since the radiated gluon 
energies are the mass splittings (of the order the Rydberg). Thus, as was pointed out in 



21| , previous attempts at using this expansion to calculate higher order corrections to quark 
mass relations p3| are on dubious ground. However, there is phenomenological evidence to 
the effect that keeping only the leading term in (|T^) may not be unreasonable for the case 
of decays (see e.g. ||12[). For the case of interest here we are concerned with self interactions 
and the OPE is a systematic expansion in Aqcd/cq- We will come back to a detailed analysis 
regarding this issue in a future publication. 

5 Mapping Onto the Chiral Lagrangian 



It was pointed on in that if a given double-dipole transition is dominated by the leading 
term in the OPE ([16D , then it is possible to map this interaction onto a chiral Lagrangian, thus 
allowing us to reduce the number of unknown parameters. Matrix elements of the double 
electric dipole operator a^-E" ■ E"", or the electric-magnetic dipole interference operator. 



^This similarity was utilized in |^. 
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agE"" ■ if" [15 , can be mapped to a chiral Lagrangian by using the fact that 



d^fn = Tj: = \^G%G^^'^ + $:(1 - 7.) m,^>.^>,. (20) 
^ 9 i 

and 

^« = ^G%G^''^ + Y: 2m,vP,75V^,. (21) 
Stt ^ . 

In the above equations, the sum is taken over the three hght quark flavors, is the axial 
current, is the dilatation current and T^i, is the stress energy tensor. 

We see that knowing the symmetry transformation properties of the Goldstone boson 
fields in the chiral Lagrangian allows us to determine how the quarkonia couple to the 
light degrees of freedom. This will allow us to disentangle the many scales which arise in 
the effective field theory. Furthermore, both operators in (^) are renormalization group 
invariants ||22||, and thus we need not worry about at which scale to choose the value for the 



strong coupling. 

Consider the case where the coupling to pions is dominated by the leading term in the 
OPE of the double-dipole approximation, namely E"^. How does this operator map into the 
chiral Lagrangian? Following (|20| ) we calculate the trace of the stress-energy tensor for the 
chiral Lagrangian which is given to order in a chiral expansion by 

df^fn = = -yTr(9^St9^'S) - 2/VTr(mS + mSt) + .... (22) 

Here the terms which are left off are down by powers of 47r/. Thus, given the analysis in 
the previous section, we now know that the lowest order coupling between the quarkonium 
state and the light degrees of freedom is given by 

9 



-'int 



-C^rlh^^^* ■ h^''^ (^d,f^ - ^(1 - 7.) m,^,^. 



(23) 



2m 

Here 7^ is the anomalous dimension for the quark field \Efj. We have neglected the term in 
the relation between E'^ ■ E"- and G^^G^^ which involves the magnetic field as it is suppressed 
by spin symmetry. 

We may now argue that the term coming form the anomalous dimensions is small com- 
pared to one based on the idea that it has arisen from integrating out shorter wavelengths 

10 



and should thus be suppressed by powers of the strong couphng evaluated at a perturbative 
scale. If we drop this term then we are left with the following coupling 

Lh'h. = ^2rlj^h^''^* ■ {fTi (9^St) + 3/VTr (mS + mSt)) + .... (24) 

If we are not willing to assume that the piece coming from the anomalous dimension is small, 
then the relative sizes of the two terms in ( p^ becomes incalculable and we must include a 
new unknown parameter into the Lagrangian. Notice that the assumption that the first term 
in the OPE of the leading twist sum dominates effectively reduces the number of possible 
terms in the Lagrangian at this order. For instance the terms 

fhlhlTrd>'Ed''i:\ K ■ KfTr{v ■ d)E{v ■ d)E^ (25) 

are higher order in the OPE. 

5.1 Derivation of mass splittings formula 



We may now use (p^ to calculate the pion mass dependence on the level splitting. The 
splitting will get a contribution from a tree level counter-term which is independent of the 
pion mass, and the leading dependence from the pion mass will come from the second term 
in (^4]). At order O(m^) there will be a one loop correction which will contribute a piece 
which is non-analytic in the pion mass as well as an unknown counter-term. We ignore 
contributions which are off diagonal in the quarkonia fields since they contribute only at 



higher orders. The mass splitting is then given by 
Am = A + Bf 



1 / TT)'^ TYl^ TYl^ 

I log + 3{m^r log + 4(m;,)^ log 



(47r/)2 47r/i2 ^ 47r/i2 47r/i 



2 



- 6C((m^)' + 



\2 



2 



+ ^ (3(m.)^ + A{mKY + {m,)') +... (26) 

The coefficients A, B, C and D are independent of the pion mass. Any pion mass depen- 
dence that one might have expected to have arisen from integrating out shorter wavelengths 
will necessarily have a well defined Taylor expansion about = and will contribute to 
a redefinition of D and other higher order chiral symmetry breaking terms. The parameter 

11 



/ is independent of the pion mass, and for consistency we must know its value up to cor- 
rections O(m^). At lowest order f = f„ and corrections to this relation can be calculated 
by computing the matrix element of the axial current between the vacuum and the one pion 
state as is usually done. There will be an unknown counter-term which, in the three flavor 
theory, can be extracted by measuring the ratio ^||25|. However, while / is independent of 



m^r, it is not independent of the number of flavors. Thus, given that all the simulations to 
date were performed in the two flavor theory, we can not determine the value of / to order 
m^. Fortunately, this is unnecessary since all the corrections will go into redefining the 
parameter D, which we are going to fit anyway. However, since we will not be fitting the 



piece which is non-analytic in the pion mass we use the relation p6 



Keeping only the leading term in the OPE, we have in addition the constraints on the 
coefficients: 

C = 1 (28) 

and 

B - (29) 

While C = 1 is model independent, we see that the numerical value of B depends on the 
model wavefunction for the onium state. 

If we reduce to the case of two families and use the relation (^) we arrive at 



Am = A, - BMJI + ^ In j-^) + (30) 
Again, when the first term in the OPE dominates, we have C2 = 1 and B2 = C2r%-^. 

6 Conclusions 

We have shown that it is possible to calculate the pion mass dependence of the levels splittings 
in quarkonia up to a few unknown constants using a chiral Lagrangian. The number of 
unknown constants which need to be extracted depends upon whether one is willing to 
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accept the approximation that the leading term in the OPE of the leading multipole result 
dominates. Within this approximation, the number of unknown parameters can be reduced 
by mapping the dilatation current in QCD onto the chiral Lagrangian. 

We have presented the a full result to order m^, which could in principle be used not only 
to calculate errors due to the pion mass, but also errors due to SU (3) breaking, should lattice 
calculations reach the point where it is possible to calculated with three light quarks. Using 
the two flavor form of our results we have calculated the errors in the extractions of ajfg{Mz) 
due to the use of unphysical values of the pion mass. Our preliminary results indicate that 
the errors in the extrapolation are of the same size as the errors quoted in ref. 0. We 
emphasize that these results are only preliminary, since we arrived at our numbers using 
data points for which our expansion in is dubious. We have found however, that under 
the assumption that the unknown counter-term at order is of the same order as the Log, 
the expansion is well behaved. Finally, it is very interesting to compare the scale f^r^ with 
the slope determined using the Monte-Carlo data. We find that both numbers are on the 
order of 10~^ MeV'-^ which we find very encouraging. This number is a factor of 10 smaller 
than the value determined using Eq. ( p9D if the model of ref. ||^ is used foifj C2. Therefore, 
at worst we have underestimated the error in the extrapolation. 
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^While we do not trust the model to predict C2 accurately, we believe the sign, which agrees with the 
one extracted from the Monte Carlo data, to be trustworthy. 



13 



References 



[1] For a status report with references see G. Alterelli, CERN-TH-95/196 , talk presented 
at Lake Louise Winter Institute, Lake Louise, Alberta 1995. 

[2] By C.T.H. Davies, K. Hornbostel, G.P. Lepage, A. Lidsey, J. Shigemitsu, J. Sloan Phys. 
Lett. B345, 42 (1995). 

C.T.H. Davies, A.J. Lidsey, P. McCallum, K. Hornbostel, G.P. Lepage, J. Shigemitsu, 
J. Sloan, |hep-lat/95 100061, Jul 1995. 8pp. Presented at International Symposium on 



Lattice Field Theory , Melbourne, Australia, 11-15 Jul 1995. 



[3] S. Aoki et. al. |hep-lat/9411058| . 3pp. Presented at LATTICE 94: 12th International 



Symposium on Lattice Field Theory, Bielefeld, Germany, 27 Sep - 1 Oct 1994. (Nuc. 
Phys. B42, Proc. Suppl. (1995) 303-305). 

M. Wingate, T. DeCrand, S Collins, U.M. Heller. COLO-HEP-355, Jan 1995. 45pp. 
Phys. Rev. D52, 307 (1995), |Ei^at/9501034 . 



[4] M. Voloshin, Nuc. Phys. B154, 365 (1979). 

[5] M. Peskin, Nuc. Phys. B156, 365 (1979). 

[6] G. Bhanot, W. Fischler and S. Rudaz, Nuc. Phys. B155 208 (1979). 

[7] P. Lepage, private communication. 

[8] J. Shigemitsu, private communication. 

K.M. Bitar, R.G. Edwards, U.M. Heller, A.D. Kennedy, P. Vranas., Talk given at LAT- 
TICE 94: 12th International Symposium on Lattice Field Theory, Bielefeld, Germany, 
27 Sep - 1 Oct 1994. 

[9] G.P. Lepage and H.B. Thacker, in Field Theory on the Lattice, Proceedings of the 
International Symposium, Seillac, France, 1987, edited by A. Billoire et. a/. [Nuc. Phys. 
B4 199, (1988) (Proc. Suppl.)] 



14 



[10] G. Brudman and J.F. Donoghue, Phys. Lett. B280, 287 (1992). 
M.B. Wise, Phys. Rev. D45, 2188 (1992). 

T.M. Yan, H.Y. Cheng, C.Y. Cheung, G.L. Lin and H.L. Yu, Phys. Rev. D46, 114^ 
(1992). 



[12 

[13: 

[14 
[15 

[16 

[17; 



[19 
[20 

[21 
[22 



[11] For a review see M. Wise |hep-ph / 9306277 



M. Voloshin and V. Zakharov, Phys. Rev. Lett. 45, 688 (1980). 

T.M. Yan, Phys. Rev. D22, 1652 (1980). 

B.L. loffe and M.A. Shifman, Phys. Lett. B95, 99 (1980). 

V.A. Novikov, M.A. Shifman, A.L Vainshtein and V.L Zhakarov, Nuc. Phys. B165, 55 
(1980). 

R. Casalbuonni, A Deandrea, N. di Bartolomeo, R. Gatto, F. Ferugho and G. Nardulh 
Phys. Lett. B302, 95 (1993), Phys. Lett. B309, 163 (1993). 



T. Mannel and R. Urech, hep-ph/95 10406 . 



G.P. Lepage, L. Magnea, C. Nakhleh, U. Magnea and K. Hornbostel, Phys. Rev. D46, 
4052 (1992). 

K. Gottfried, Phys. Rev. Lett. 40, 538 (1978). 

Weak Interactions and Modern Particle Theory, H. Georgi, Bejamin/Cummings Menlo 
Park Ca. (1984). 

M. Luty and R. Sundrum, Phys. Lett. B312, 205 (1993). 

A. Zee, Phys. Rev. Lett. 29, 1198 (1972). 

J. Lowenstein and B. Schoroer, Phys. Rev. D7, 1929 (1973). 

S. Adler, J.C. Colhns and A. Duncan, Phys. Rev. D15, 1712 (1977). 

B. Grinstein and L. Randall, Phys. Lett. B217 335, (1989). 



15 



[23] J.F. Donoghue and D. Wyler, Phys. Rev. D45, 892 (1992). 

[24] Y.P. Kuang and T.M. Yan, Phys. Rev. D24, 2874 (1981). 

[25] J. Gasser and H. Leutwyler, Nuc. Phys. B250, 465 (1985). 

[26] P. Langacker and H. Pagels, Phys. Rev. D8, 4595 (1973). 



16 



